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1. Introduction 

It is well known that the BFKL Pomeron calculus |l], ^ gives the simplest and the most elegant approach to 
the high energy amplitude in QCD. This approach is based on the BFKL Pomeron |^] and the interactions 
between the BFKL Pomerons [Q, ^, ^, Q , which are taken into account in the spirit of the Gribov Reggeon 
Calculus P] . It can be written in the form of the functional integral (see Ref . B ) and formulated in terms 



of the generating functional (see Refs. |^, IC, 11, 12, |T^). However, in spite of the fact that we have learned 
a lot about the properties of the high energy asymptotic behaviour at high energies, even the simplest case 
of the BFKL Pomeron in zero transverse dimension (a toy model , see Refs. ||9|, ^, |T^) has not been solved^ 
In this paper we develop the analytical approach to this problem and find the solution. We hope that this 
solution will help us to develop a technique for the general approach in which the amplitude depends on 
the size of interacting dipoles as well as on their impact parameters. 

The simple toy model has some advantages which makes it a good laboratory for seeking the solutions 
that include the Pomeron loops. First of all, the mean field approximation ||l|, |2|, |l^, 17] for this model has 



a simple analytical solution |T^, ^] while it is a difficult problem in the general case and we have only few 
analytical approaches |18, ^], therefore, we have to rely very often (too often) on the numerical solutions 



20[ . The JIMWLK approach ||31|] coincides with the Balitsky-Kovchegov ||r^ one in this model making 
our life simpler in the mean field approach. 



The first corrections due to the Pomeron loops which were suggested in Ref. [ESl M can be easily 



summed in this model |3C] giving a guide for the general impact of the Pomeron loops on the high energy 
asymptotic behaviour of the scattering amplitude. 

We are aware that we can loose the important property of the QCD approach (see Refs. [32, 53, |3, 



35, pm) but we believe that the analytical solution to the simple model is a necessary step in finding more 
general and certainly more complicated solutions. We are certain that methods developed for solving this 
model, will be useful for searching for a general solution of this challenging and difficult problem. 

The paper is organized as follows. In the next section we discuss the general properties of the equation 
that we need to solve. Being a diffusion equation this equation could be rewritten as the Sturm-Liouville 
equation which has a discrete spectrum. In section 3 we develop the semi-classical approach which results 
in finding the analytical solution to the equation. We show that this solution approaches the so called 
asymptotic solution at high energies. In this section we discuss the analytical solution at u ^ 1 which 
gives the transparent example of the main properties of the general solution. 



^Sortie progress has been made in Refs. |l3, u4, UM and we will comment on these efforts below. 



In section 4 we find the general analytical solution to the problem. We state that the prolate spheroidal 
wave functions is the complete set of the eigenfunctions of the master equation and build the Green function 
as well as solution to the master equation that satisfy the initial and boundary conditions that has been 
discussed in section 2. 

In conclusions we discuss our results and check the approximate approaches, developed for this kind 
of equations, against the exact solution. 

2. The equation and its general properties 

2.1 The equation and initial and boundary conditions. 

In the toy model the equation for the generating functional^ has a simple form 0, 12, O] 



_= „(!-„) (^-._ + _j (2.1) 

where 

y = T{2^l)Y = {T{1^2)/k)Y (2.2) 

with vertex r(2 — > 1) which describes a merging of two dipoles into one dipole. k is the large parameter of 
our problem and it is equal to 



2 • r(l ^2) _ 1 

?5 



- - -^n-0 — T^ - — » 1 (2.3) 

i (z ^ Ij ai 

where as is QCD coupling and r(l — > 2) = {Nc/'K)as = as is the vertex for the decay of the dipole 
into two dipoles. The order of these vertices r(l — > 2) and r(2 — > 1) are chosen from the general equation 
for the generating functional in QCD (see Ref. ||l^ for example). 

Eq. (|2.l| ) is a typical diffusion equation with the diffusion coefficient that depends on variable u. We 
need to add an initial and boundary condition to solve such kind of equation. Generally, this equation 
describes the interaction between Pomerons with intercept Ap = r(l ^ 2) = as- The Pomeron diagrams 
that describes Eq. ( |2.lD are shown in Fig. |l[ 

In this paper we are going to calculate the sum of enhanced diagrams (see Fig. |l|) and for them we 
have the following initial condition ^ 

Z{Y = Q;u) = u (2.4) 

Since the physical meaning of Z is 

oo 

Z{Y-u) = Y, Pn{Y)u^ (2.5) 

n=0 



■^Actually for the toy model, in which there is no dependence on the sizes of interacting dipoles, the generating functional 
degenerates to the generating function. 



In Eq. (2.4) as well everywhere below we use the notation Y for y. Hopefully it will not lead to troubles in understanding. 



r{i->2)Y 



r(l -> 2) 

-r(2->i) 




Pomeron and vertices of Pomeron interactions 




enhanced diagrams 



mean field approximation 



Figure 1: The Pomeron diagrams that are summed by Eq. (2J). Sohd Hnes denote the Pomerons. The first three 
diagrams show the Pomeron and vertices of Pomeron interaction. The next five diagrams are the examples of so 
called enhanced diagrams which we are going to sum in this paper. The last two diagrams give the example of the 
Pomeron diagrams in the mean field approximation. 



where Pn(X) is the probability to find n-dipoles (or/and n-Pomerons) at rapidity Y. Therefore, Eq. (2.4) 
means that at low energies we have only one Pomeron. 



The master equation (see Eq. (2T)) describes Markov's chain for the probabilities P]\[{Y), namely 
dPn{Y) 



dV 



-r(l ^ 2) n Pn{Y) + r(l ^ 2) (n - 1) Pn-i{Y) 

- r(2 ^ 1)^^^ PniY) + r(2 ^ 1)!^^^ p„,,(y) 



Eq. ( |2.6| ) has a simple structure: for every process of dipole splitting or merging we see two terms. The 
first one with the negative sign describes a decrease of probability Pn due to the process of splitting or 
merging of dipoles. The second term with positive sign is responsible for the increase of the probability 
due to the same processes of dipole interactions. 

The boundary conditions are very simple 

Z[Y-u = \) = J2 Pn{Y) = 1 (2.7) 

n=0 

which directly follows from Eq. (|2.5| ) and express the conservation of probabilities, saying that the total 
probability is equal to unity at any value of rapidity. 

Using the Laplace transform 

r-a+ioo ^^ 



ZiY;u) 



we can reduce Eq. {^^ to the following equation 



a— too 
4 



27ri 



Z{ll!;u) e' 



jY 



UJ 



u{l — u 



Z{u)]u) = -kZ'^{u)]u) + Z" {u:]u) 



(2. 



(2i 



''in principle, this equation contains an inhomogeneous term (initial condition) Z(Y — Q,u) , iiowever, we will include the 
correct initial condition latter, after solving this homogeneous equation. 



Introducing new functions 



siu) = — re"''"; and p{u) = e''"'' ; (2.10) 

u(l — u) 



we can reduce Eq. ( p.9| ) to the typical Sturm-Liouville equation [21, |2^ 

— s{u)uj Z{uj,u) + — {p{u) Zu{io , u)) = (2-11) 

2.2 General properties of the Sturm-Liouville problem 

It is very instructive to remind ourselves of the general properties of the Sturm-Liouville problem which 
has been studied by mathematicians for a long time [^, 22]. 



1. Eq. ( 2.11| ) has the infinite set of the eigenvalues iOn = — A^. All A„ are real and can be ordered so 



that Ai < A2 < A3 with A^ — > +00. Therefore, there can exist only a finite number of negative A^. 



But for our particular case, (for Eq. ( 2.11 )) there are no negative eigenvalues An and Ai = is the 



least eigenvalue, to which there corresponds the eigenfunction Zq{u) = Constant; 

2. The multiplicity of each eigenvalue is equal to 1; 

3. The eigenvalues are determined up to a constant multiplier. 

4. Each eigenfunction Zn{u) has exactly n — 1 zeroes in the interval (0, 1); 

5. Eigenfunctions Zn{u) and Zm{u) are orthogonal with weight s{u)^ namely, 

/ dus{u) Zn{u) Zm{u) = for n ^ m (2-12) 

Jo 

6. An arbitrary function F{u), that has a continuous derivative and satisfies the boundary conditions 
of the Sturm-Liouville problem (in other words, a function that we, as physicists, want to find), can 
be expanded into absolute and uniformly convergent series in eigenfunctions: 

oo 1 /"^ 

F{u) = Y^PnZniu); where F„ = --—r / du' F{u') Zn{u') (2.13) 

1 ll-^nll Ju 



n=l 

where 



\Znf = / dus{u)Z^{u) (2.14) 

Jo 



7. The following asymptotic relation holds for large eigenvalues as n — > oo: 



A. = ^ + 0(1); where A = £ J ^ ^^ ; (2-15) 



8. In our case (see Eq. ( 2.10| )) 



f du' 

A = / , = TT, 2.16 

Jo ^/u'{l - u') 

and, therefore, at large n 

Xn = r? (2.17) 

In Appendix E we compare this asymptotic spectrum with the spectrum of the exact solution as well 
as with the perturbation procedure to calculate it, given in the next section. 



9. In out case we can use Eq. (2.17|) and Eq. (|2.15) for n ^ 1 to estimate not only the values of \n 



but also to find the corresponding eigenfunctions. These eigenfunctions have the following form: 



,^ „ - ,^m(1-u) 62" sin(nG(u)) + O - (2.18) 

\Zn\\ V^ / \'^) 

where 

<d{u) = - C ff = -f - arcsin(l-2u) (2.19) 

Ju ^Ju'(\-U') ^ 



Introducing a more general equation ^ than Eq. ( ^.11 ), namely. 



where ^iY, u) is a given function. We can find a general solution to this equation with the initial condition: 

Z(y = 0,n) = f{u) (2.21) 

and with the boundary conditions: 

Z (y, u = 0) = 0; (2.22) 

Z{Y,u = l) = 1; (2.23) 



The form of this solution is (see section 1.8.9 in Ref. [p2| ) 

Z(Y,u) = j dV j d5<i>(y',f) G{Y-Y':u.i) + j d(4li)f{()G{Y;u,() 

Jo 9^ 

+ p{u = 0) r dV Z{Y,u = 0) ^^-^--p^^l^^, (2.24) 

Jo (^^ 



'More useful information on the general Sturm-Liouville is given in Appendix A and in Refs.M, 02 



where the Green function of this equation is equal to 

G iY;u,0 = f: ^" ^r^ ^; ^^^ e-^"^ (2.25) 

n=l 

where Z^ are the eigenfunction of Eq. ( p. 11 ) with the boundary conditions 



Z„(u = 0) = 0; Zn{u = l) = 0; (2.26) 

One can see that we know a lot about the properties of the master equation. 

2.3 The iterative algorithm for a search of the solution 

The general formula, given in the previous section, leads to a definite algorithm for numerical calculations. 
In the first stage of this algorithm we can use the eigenvalues given by Eq. ( |2.17| ) , and the eigenfunctions 
from Eq. ( p. 18 ). The Green function has the form 



Goo{Y;u,^) = - Vsine(n) sinG(0 gt ("+«) V exp(in{e(u) - 9(0}) e""'^ (2.27) 

n=— oo 

The solution can be found from Eq. ( |2.24| ) and it has the form 

Zoo {Y, n (6) ) = r clY' f clQ' ^(^°^Q^+ ^) g— O) q^ {Y - ¥'; n (9) , ^ (9')) - (2.28) 

Jo J-K sm9' 



di{Q' 



dY Tr-TTTT. |e'=o 



In the next stage of the algorithm we are searching for the solution in the form 

Z{Y,u) = AZ(y,u) + Zoo{Y,u) (2.29) 

and for AZ (a;,u) we obtain the following equation 

- s(u)ojAZ{oj,u) = L{AZ{uj,u)) + L(Zoo(u;,'u)) (2.30) 

where we denote by L the Sturm-Liouville operator, namely L = -^ ip('^)du)' 

Eq. ( p. 30 ) is prepared for the iterative procedure of finding the solution using Eq. ( |2.28| ) as the first 



iteration. To calculate the second iteration we substitute Z^o (Y, u) into Eq. ( |2.3C| ) which has the form of 
Eq. ( |2.2C| ) with ^{Y, u) = L {Z^o {Y, u)). Using Eq. ( |2.24| ) we can find AZ {Y, u). Repeating this procedure 



several times we will obtain the numerical solution to the problem. However, the procedure would be better 
converged if we could determine the eigenvalues of the problem at small values of n. Below we suggest two 
approximations that allows us to approach small values of n. 



The algorithm described above has a very simple meaning if we re- write our equation as a Schroedinger- 
type equation. It is easy to see that the search for Zqo {Y, u) of Eq. ( |2.28| ) can be simplified if we introduce 
a new function T 

Z{Y,u) = T{Y,e) et« (2.31) 



For this function we can rewrite Eq. ( |2.9D in the very convenient form for the function T (Y, G) intro- 
ducing a new variable (see Eq. ( p.l9D ) instead of u. 

" ^(-,0) = -^sinGr(.,e) + 4,-^^^^ (2.32) 



sine ' ' 4 ' ' ' dGsine dQ 

which is a generalized Sturm - Liouville equation of the form 

,(e)im^^A„e)^^I^-„e)Ta,e) (2.33) 



with 



k2 



s(G) = l/sin0; p(0) = l/sinG; g(0) = — sinG (2.34) 

Eq. ( 2.33| ) allows us to find corrections of the order of 1/n^ to the spectrum of the Sturm - Liouville 
operator. Indeed, A„, is equal to (see section 1.8.9 in Ref. [p2|) 



= n + —Q{0,-tt) + o(\] (2.35) 

nn yn-^ J 

where 

1 r® «-2 

Q(e,Q') = - / dQ" q{Q") = — {cosG' - cosG} (2.36) 

The estimate for the eigenfunctions looks as follows 

r„(0) = sin(n0) - -^ {2G + ^(1 - cosG)} cos(nG) + O ( ^ ) (2.37) 

47rn \n'' J 

Eq. ( 2.32| ) can be re-written as a Schroedinger-type equation: 



(dQ) 



P)\2 



{E-U{Q))^ = 0. (2.38) 



Indeed, Eq. (|2.32D has the equivalent form 



^2 



ojT = -— sin^GT - cot G Te + Tee (2.39) 

Introducing (see Appendix C for general way of reducing the Fokker-Planck equation to the evolution 
described by hermitian Hamiltonian) 

T = (/.(0)^(G) with (/.(G) = VsinG (2.40) 



we obtain for ^ the Schroedinger equation with 



E 



-LO 



and 



(2.41) 



c/(e) 



120 



100 




— sin^e + - (2 + 3cot2e) ; 

4 4 

In Fig. ^ we plot this potential. One can see that 
for large k it has a typical two minima form with the 
maximum at = tt/2 while at — > and — > tt 
f/(0) increases approaching infinity. 

Therefore, for large energy excitation we can re- 
place this potential by the rectangular-well potential 
with the wave functions given by Eq. ( 2.18| ) with the 
spectrum of Eq. ( 2.15| ). Indeed, Fig. @ shows that U {Q) — 



sua. 



(0) 



(2 + 3cot^ 0) denoted as dashed line 



, can be replaced by the potential which is equal to 
zero for < < 27r and which is infinitively large 
outside of this region. Such potential allows us to re- 
duce Eq. ( 2.38| ) to the Schroedinger equation with the 
simplified potential, namely, 

' ~ -^^Q)]^ = 0; with ^ (0 



Figure 2: The potential energy U (6) of Eq. (|2.4lD 
versus 8 at k=16, . The dashed line shows the po- 
tential U (Q) — {k'^/4:) sin'^ 9 which corresponds to 



(d0) 



P)\2 



+ E 



K 

T 



sm 



0) 



^ 



the set of the eigenfunctions given by Eq. ( 2.18 ) 



(2.42) 
The boundary conditions in Eq. ( 2.42| ) reflects the fact 
that potential is very large outside the [0, 2tt] interval 
in 0. 

The fact that we have two minima potential leads to the asymptotic behaviour of the amplitude which 

for w = (see Refs.p, P, 



is determined by the asymptotic solution of Eq. 



I), namely, 



Z, 



asymp 



{Y = CO, u) 



1 



(2.43) 



which satisfy both boundary conditions: Zasympiu = 0) = and Zasymp{u = 1) = 1 . 

In section 4 we will give the exact analytical solution to the master equation and we postpone until this 
section our discussion of this Schroedinger- type of the equation. However, this equation (see Eq. ( 2.38[ )) 
allows us to use the usual perturbative approach for the numerical solutions. For example we can represent 
the potential C/(0) of Eq. (|2.41D in the form 



C/(0) = C/o(0) + y(0) 



(2.44) 



and solve Eq. ( ^.38 ) with Uq exactly. Considering V (0) as being small we can develop the usual pertur- 
bation approach in quantum mechanics for a numerical solution of the equation. The value of k is large 



in our approach and we need to take such a perturbation procedure with a great precaution. In appendix 
E we give the spectrum of A„ calculated in the first order of preturbation approach and it terms out the 
agreement with the exact spectrum is very impressive. 

However for large values of k the more reasonable way to find a good first approximation for searching 
Zoo is to use the semi-classical approach which we consider below. 

3. Analytical solution to the Sturm-Liouville problem 



3.1 Semi-classical approach 

In our master equation (see Eq. (|2.9| )) we have a natural large parameter: k (see Eq. (|2.3D ). We wish to 
use this parameter to develop a semi-classical approach to the master equation. We assume that 



Z{u>] u) 



("^'"^ where 



<C 



Substituting Eq. (3.1) into Eq. (2.£) one can reduce this equation to the form 



n(l — u) 



'l^<Pu + ('/'«)' 



(3.1) 



(3.2) 



Solving Eq. (3.2) we obtain 



^i 



du 



K I / 4a; 



and 

(fr^ {u;u) = - I 4>^,{lo;u') du' 

J u 

The general solution to Eq. ( ^.1| ) has the form 



Z{Y-u) 



a+ioo 



duj 



i^^-\uo)e^ (^'") + $(+)(w)e'^''('^'")} 



(3.3) 



(3.4) 



(3.5) 



where functions $'- •'(w) and ^^'^'{uj) should be found from initial and boundary conditions (see Eq. (2.4) 
and Eq. <^). We recall that we use Y ioi y = T{2 -^ 1) Y. 



The contour of integration (contour C in Fig. ^) over u> is situated to the right of the singularities as it 
is shown in Fig. ^. As we have discussed the general structure of the singularities looks as follows. We have 
no singularities for positive uj = u)^ > 0. It looks strange since one can see from Eq. ( |3.3D that we have 
a restricted number of singularities for positive uj = lo^ > 0. However, since the contour of integration 
C in Fig. m is placed between loq = and to = ujf, we can close it to the left semi-plane (contour Ci in 
Fig. ^ and, therefore, the singularities cj+ does not contribute to the integral of Eq. (p.Sp. 
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•"© 



COn 



co: 



co^ co^ 



COrf - ?^n <0 



C0| (JO2 03+ 



The analytical function for Eq. (|3.4| ) is writ- 
ten in the Appendix B while here we will discuss 
this function in two different limits: ■> f^ — r <C 
1 and -^r^^^ » 1- 

For -5— 7T — r ^ 1 we have for function S^ (cj; xl) 
the following expression 



max 






UJ 



k{u — 1) + — In 



l-u 



LO 



In 



u 



(3.6) 



(3.7) 



It is easy to see that (j) {uj; u) of Eq. (p.tjj ) cor- 
Figure 3: The contours of integration over uj in Eq. (|!|). responds to the solution which survives at k ^ 
All singularities are located along real negative real axis ^ ^^-^^ describes the case of the mean field ap- 

proach. Analyzing the mean field solution (see 
Ref. |]lO[ , for example) one can see that all positive singularities u;"*" in Fig. y) are the singularities of $'~) 
while $'+) has the singularities at the negative values of to. Therefore, the only singularities that we need 
to take into account are those that are described by Eq. ( |3.6| ). 

For u; ~ K^ we need to expand the exact formulae for (j) {lo,u) in the kinematic region. For 

K^uii^u) ^ 1- In this region (j}^{uj,u) are equal to 



<P^{^,u) 



(3.8) 



;{u-l) ± iVl^l Q{u) 



'-{u-1) ± i^n Q{u) 



where is given by Eq. ( 2.19| ) and w„ = —^n = — Mn- 

Therefore, in both branches ^'^ we have singularities at negative w and from estimates of Eq. ( |2.17 ) 
we know that A„ = n. Using simple relations between variables u and G 



u = - (cos + 1) 



we can rewrite exp ((/)^(u;,u)) in the form 



exp 



'-{U)n = -K = -f4i^u)) 



where we introduce new variable 



c 



K 

exp I - 



ie 



-2+C+ ^ 



"i/in 



(3.9) 



(3.10) 



(3.11) 



It is obvious from Eq. ( |3.8| ) that C" is the eigenfunction of our master equation at large values of u 
((J > k). However, the eigenfunctions that satisfy the boundary conditions of Eq. ( 2.26| ) are 



Zni@) = sin(nG) 



(3.12) 



11 



These functions form the complete set of functions in the interval < n < 1. In other words, the 
approximation of Eq. ( p.l7| ) turns out to be exact in our semi-classical approach. Having this fact in mind, 
we can find function <l>(^-'(u;) in Eq. ( |3.5D from the initial condition of Eq. ( |2.4| ). 

Functions C" are eigenfunctions of Eq. (2.32) in which we neglect the first term in the r.h.s. and the 
second one is modified assuming that 

d 1 dT(io,e) . , ,, 1 ^^^-(^^e) 

IS replaced by 



desinG d& "^ ' sine dQ'^ 

In this approximation Eq. (p.32|) reduces to 



d'T{u:,Q) 
u;T{uj,Q) = —2 (3.13) 



dQ 



It is easy to see that Eq. ( 3.12| ) gives the eigenfunctions of this equation. 



Using the solution of Eq. ( p. 13 ) we can calculate the Green function for our master equation (see 



Eq. ( |2.9D ) that this Green function has a very simple form, namely 

G^ {Y;e,Q') = iet+t[=°^0 + ™«e'l ^ sin(ne) sin(ne') e""'^ (3.14) 

n=l 
and the initial conditions have the form 



r(y = 0,e) = -(l + cosG) exp('-j(l + cosG)) (3.15) 

Using this Green function we obtain the semi-classical solution in the form 
Z^^ = (3.16) 

- Yl S e-"'^et" sin(nG) / (iQ'sin(nQ') - (1 + cosQ') exp (- ^ (1 + cosG')) [ 

n=-oo '^ '' ~'" ' 

The integral over G' could be taken using the generating function for the modified Bessel function of 
the first kind (see formula 9.6.33 in Ref.||2f 



exp(^cosG) = Yl ^"(l) e^^"" (3-17) 

n=— oo 

and the expression for simple integral 

k"^ — n? 



dQ' cos(nG') sin(fcQ') = - ^ J^ i '- (3.18) 



From Eq. ( |3.17| ) and Eq. ( 0.18 ) one can derive that 



C{n,K) = f dG' sin(nG')^(l + cosG') exp(-j(l + cosG')) =^exp(-K/4)x (3.19) 

J—n 2 V 4 / 2 
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Y, (/™+i(-^/4) + /„_i(-Ac/4) + 2/„(-k/4)) ^ ^,^ J, ^ + (/i(-k/4) + 2/o(-^=/4)) ^ ^ 



^m=l 



■n? — m? 



n 



Substituting Eq. ( 3.19 ) into Eq. ( |3.16| ) we obtain 

Z^ = iet("-i) |C(0,A:) + 2^ {e-"'^et" sin(ne) C7(n, k)} 



TT 



(3.20) 



n=l 



However, this expression for Z^ has two disadvantages: (i) we approximate uJn = —v? even at 
small values of n where we know that our spectrum is different (see Eq. (|3.6|) and Eq. (^.7|)); and (ii) 
Z^{Y,u= 1) / 1 (actually Z^{Y,u = 1) = 0). 

The set of eigenfunction at small values of n is clear from Eq. (|3.6[) and Eq. (|3.7D , namely, 



Zn{u) 



for u > 12 



for u < 2 



pK (m— 1) ^n 



r" = exp f K j^ I r"' where r 



M-l 



(3.21) 



r" where r 



u-l 



The fact that we have two different sets of the eigenfunctions for u > 1/2 and u < 1/2, directly comes 
from the contour integral of Eq. ( |3.5D and the expressions of Eq. (|3.6| ) and Eq. ( |3.7D for functions 0^ and 
4>'^ . Indeed, for u > 1/2 ln(u/(l — u)) is positive and we can close the contour in Eq. ( p.5| ) to the left 
semi-plane for 0"^. For u < 1/2 this log is negative and we have to use 4>~ for negative values of tu. Recall, 
that we know from the general properties of the Sturm-Liouville equation that only negative uj contribute. 

The value of iOn = —An for the eigenfunctions of Eq. ( [j.21 ) is equal to 



An 



K n 



Coefficients C{n, k) we can find using the following series 



u exp (— k(m — 1)) 



1-r 



exp — K 



n=0 



(3.22) 



(3.23) 



where Ln{—K) is the Laguerre polynomial (see formula 8.975 in Ref . pq] ) . 

From Eq. ( |3.23| ) one can see that the initial condition Z{Y = 0, u) = u can be re- written in the form 



of 



Y, Cin,K)Zn{u) q(u>\) 

n=0 ^ '^^ 



(3.24) 



withC(n,K) = Ln{—K). 

We suggest to build the Green function considering for all A„ < k? the set of the eigenfunction given 
by Eq. (|3.21| ) while for A^ > k^ we use the set of Eq. (|3.12|) . 
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It means that our Green function has the form 



n=\K.]+l 

n kY 



G^^(y;7x,0 = exp(K(^ + e-2)) ^ T{u)T{i)e 



n=0 

oo 



+ exp('^(^i + ^-2)) Y^ sm (ne(n)) sin (nG(e)) e""'^ (3.25) 

n=[K]+2 



which leads to the generating function in the form 



n=[K]+l 



( , ,1,— ^rvj- 

u > -\ exp(K(u-l)) Y^ 
^ n=0 

n=[K]+l o<j 

+ eU<-j Y r(n)"e-""^ + exp(|(u-l)) ^ sin (nG(u)) C(n, k) e 



t(u)"L„(-k) e""'^^ (3.26) 

n=0 

n=[K]+l 



-n^y 



n=0 n=[K]+2 

This solution has two problems which cannot be solved analytically: a dependence of the separation 



parameter between large and small n which was taken in Eq. ( 3.26 ) to be equal to [k] + 1 and rather com- 



plicated expression for C{n,K) (see Eq. ( 3.19 )) which does not allow us to make all calculation analytically. 



We investigated both problems computing Eq. ( 3.26 ) numerically for different separation parameters. It 



turns out that the answer is very insensitive to the exact value of the separation parameter if only we take 
it around k. Actually, the value of this parameter we can find comparing exact semi-classical functions cj)^ 
with our approximate ones ( see Appendix B, where we discuss the exact functions). 

The behaviour of the scattering amplitude is illustrated in Fig. ^ as a function of rapidity Y. In Fig. ^ 
one can see the u dependence of our solution at large values of rapidity Y. This behaviour confirms the 
expectation that has been discussed in Ref. |Q, that at high energy our master equation predicts the gray 
disc behaviour for the scattering amplitude. 

We check the accuracy of the semi-classical solution by calculating the following function: 

^(^'^) = ^ mmimal separate term of the numerator ^''''^ 

As we have expected the semi-classical approach is rather bad at small values of y < 0.05 leading to 
A(Y,u) ~ 40% for u ~ 0.5. However, for large values of y (y > 0.1 the accuracy is better than 10 % 
for all values of u. It should be stressed that for u > 0.8 in the entire kinematic region of Y the accuracy is 
not worse than 10%. All estimates were done for k = 16. This fact is very encouraging since the scattering 
amplitude N can be calculated in the following way |17, 0] 



iV(y,7(yo)) = 1 - Z(y,n = l-7(yo)) (3.28) 
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where 7(io) is the scattering ampUtude at low energy (y=0). As it is discussed in Refs. |17, 11] in 



Eq. ( 3.2^ ) we assume that n dipoles interact with the target independently. It means that the scattering 
amplitude of n dipoles with the target 7n(Yo) = 7"(^o)- ^ In our QCD motivated model 7 ss 1/k <^ 1. 
Therefore, we are interested only in the values of u that are close to unity (u ^ 1), where the semi-classical 
approach works quite well. 
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Figure 4: The dependence of the scattering ampli- Figure 5: The asymptotic dependence of the scat- 
tude versus rapidity at fixed k — 16. Recall that here tering amplitude (see Eq. ( 3.26| ) ) versus u at fixed 
we denote Y = y (see Eq. (2.2). k = 16. 



However, in the next section we develop a better analytical approach for such u. 
3.2 Exact solution at u — ^ 1 



To check that the eigenvalues in the region of small lo is determined by Eq. ( 3.22| ) we develop a different 
analytical approach. In the vicinity of u = 1 we can replace Eq. (|2.lD by the following equation: 

(3.29) 



dZ , f dZ 

dy = ('-"K""^^ + "(5n)2 



(3.30) 



For Z {lo,u) Eq. ( |3.29| ) reduces to the form 

ujZ{ll!,u) = —k{1 — u)Zu{uj,u)+ u{1 — u) Zuu {^ , u) 

The solution to Eq. ( 3.3C| ) is the Gauss's hyper geometric function 

Z{uj,u ) = u^+'' 2Fi{a,(3,j,u) (3.31) 

^This contribution for low energy interaction of n dipoles with target is proven in Ref. [M. He has also been used in the 
mean field approximation IitI. Nevertheless, it looks as an additional assumption for dipole-dipole scattering. We thank 
Al. Mueller and E.Iancu for very fruitful discussions on this subject. 
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(see formula 9.151 in Ref. [^) with the following parameters 

1 .. . 1 



a = - (1 + k) + - ^(1 + k)2 - 4w ; 
1 . . 1 



/3 = 2 (1 + '^) - 2 \/(l + ^? - 4u; ; 

7 = K + 2 ; (3.32) 

The generating function that satisfies the boundary condition Z (Y, u = 1) = 1 is equal to 

Z{io,u) = u^ ——. — 7 (3.33) 



Since (see formula 9.122 in Ref. H) 

r(7) r(i) 

2Fi(a,/?,7,^ = l) = jr. , p, ' ^, (3-34) 

r(7-Q) r(7-/3) 

The eigenvalues stem from the poles of r(7 — a) or r(7 — /3), namely, from the equation 

7 — a = — n ; which gives u;^ = —^n = —n — n(l + K) where n = 1,2,3,... (3.35) 

This equation leads to a smooth transition between large and small values of n supporting our semi-classical 
formula. Substituting Eq. ( |3.33| ) in Eq. ( p.5| ) one can find the eigenfunction of our equation. They have 
the following form 

^«W = Yhvl "'^" ^n''''^'(2t^ - 1) (3-36) 

where Pn ' '^ are the Jacobi polynomials (see sections 8.960 - 8.967 in Ref. |2^]). 



The normalization of function Z„ in Eq. ( 3.36 ) is obvious since 

\\Znf = (3.37) 

C du ^ (p-'^'+^(2u - D? = r(n) r(n + 2 + K) ^ (^+1 + ^) 

Jq l-n ^ " ^ '^ n!(l + K + 2n)r(l + K + n) n(l + K + 2n) 

It is easy to check that eigenfunctions of Eq. (|3.36| ) satisfy the following boundary conditions; 

Zn{u = 0) = ; Zniu = 1) = (3.38) 



We found that the easiest way to solve the master equation (see Eq. (|3.29D ) is to introduce Z{Y,u) 



u + Z{Y,u). The equation for Z can be reduced to the Sturm-Liouville form, namely, 

. ^dZ d , ^ dZ ^, ^ , ^ 



where 



s(u) = — -u '^ ; p{u) = u '^ ; $(n) = -ku '^ (3.40) 

u[l — u) 
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The solution to Eq. ( 3.391 ) has the form 
Z{Y,u) = 
where 



r dY' [ d^G{Y-Y';u,^) ^0 
Jo Jo 

= Z^ iT^Tlv?. exp(-A„y) 



n=l 



\7 l|2 



(3.41) 
(3.42) 



where A„ are given by Eq. ( |3.35| ) and Z„ is the set of the eigenfunction shown in Eq. ( |3.36| ) . 

In Fig. ^ we plotted the energy behaviour of the scattering amplitude for different values of u that are 
close to unity. One can see that this solution gives the amplitude which is quite different from the semi- 



classical solution. It sounds controversial, but actually it is not, since the solution of Eq. ( 3.41 ) contains 
the integral over ^ from to 1 while our functions Z„ of Eq. ( 3.36| ) we know only at n — > 1. 

It turns out that we can find a very simple solution at small values of u (u — > 0). Introducing a new 
variable 7 = 1 — -u we can rewrite our master equation in the form at 7 — > 1 

, f dZ d^z 

(1 -7) ^ ^r +7 



dZ 

or in (^-representation it looks as 

ujZ{uJ,-f) = (1-7) {k Z^{uj,-f) + jZ^^{uj,j)) 

The solution to this equation has the form 

Z{uj,-f) = 2Fi{a,(3,K,-f) = 2Fi{a,(3,K,l-u) 
where 

a = ^ (^-1) + V(l-'^)^ - 4c 
P = 



(3.43) 



(3.44) 



(3.45) 



[UJ 



\{{^-i) - Vo^W 



Auj 







The condition that gives us the spectrum of to is 

Z{u!,l - u = 1) = 2Fi{a,P,K,l) 

The equation for the eigenvalues has the form 

which gives w„ = — A^ = ?^ + n[K + 1) + 1 where n = 0, 1, 2, 



r(^) r(i) 

T{K-a) T{k-P) 



a 



n 



(3.46) 



(3.47) 



(3.48) 



One can see that this spectrum practically coincides with Eq. ( p. 35 ). Therefore, we can assume that 
Eq. ( 3.36 ) or Eq. ( 3.35| ) determines the spectrum of our problem in the entire region of u. This fact we 
will use below to develop the numerical approach to the master equation. 
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3.3 Approximate method 



The semi-classical solution we can compare with a simpler approximation suggested in Ref. |12]. This 
approach consists of two steps. The first one is to find the asymptotic solution, namely, the solution to the 
following equation 



^(Pu + ((PuY 
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Figure 6: The scattering amplitude as a function of rapidity 
at different values of w ^ f . Recall that here we denote by 
Y = V (se e Eq. ([2.2[)). k, is equal to f6 in this picture. 



Eq. §M) .^ 

The solution has the form 

Z{u;Y) 



ue 



f^y 







This solution is 



(3.49) 



'asymp 



(n) 



K (n— f) which 
satisfies the normalization (j){u = f) = that 
follows from Z{u = \,Y) = 1. The next step 
is to find a solution in the form i;^(u; y) = 4>asymp{u) + 
A(/)(m; y) assuming that A0(u; y) ^ (pasympiu). 
The equation for A0(n; Y) is the typical Liouville 
equation 



dA(f){u;Y) 



K u{l — u) 



dy 

with the initial condition 



dA(t>{u;Y) 
du 



(3.50) 



A(/.(u;y = 0) = \n{u) -k{u-1) (3.51) 

The solution is an arbitrary function of the fol- 
lowing type 

1 



A(/)(u;y = 0) 



H[Y + 



In 



1-u 



(3.52) 
which can be found from the initial condition of 



1 — u + ue'^y 



exp [ Ku{l 



1 



o^^y 



1 — u + ue'^y 



(3.53) 



This solution corresponds to the general solution, given by Eq. ( 3.26| ), but only one solution (l)'^{u>,u) 
is taken into account. In spite of the fact that the full semi-classical solution looks quite different from the 
simple expression of Eq. ( |3.53| ) the both solution are in less than 1% agreement at least for 3^ > 0.3. 



4. The exact solution to the problem 



A master equation of Eq. ( |2.9D we can re- write in more convenient form, introducing a new function Q, 
namely ^, 

Z {uj,u) = e^'^Q{v,Lo) u{l-u) 



''One can see in Appendix C that Eq. ([4.l|) is a particular choice of the general transformation. 
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l-f2 



g{v,u;) et^^-'^) 



4 
where v = I — 2u (4.1) 



For the function Q we have 

ug{oj,v) = -^(1 - i,')giu;,v) + ((1 - v')giu;,v))^^ (4.2) 

The equation for g can be re-written in the form 

(l-t;2)g„„ -4r;g„ + j-u; - 2 - ^(l-t;2)| g = (4.3) 



For function Sn,m = (1 ~ ''-' ) ^ g{uj,v) Eq. (O) reduces to the equation for the prolate spheroidal 



wave function with the fixed order parameter m = 1 and arbitrary degree parameter n which is an integer, 
Sn,m{c,v) (see Refs. ^, |3|), namely, 

(4.4) 
av \ uv / \ 1 — v~ / 

where 



d 


( ^, dS ic,v)\ ^ / ^2^2 
V dv J \ 


m2 \ 


5'„,m(c,i;) = 


= 


dv 


l-t;2y 




Wn = -3^ - 2 + m(m + l) - A™ 


2 
~ ~ 16 








n ft, 

m = 1 ; c = . 









(4.5) 
16 ^ ' 

To verify Eq. ( ^ ) we can substitute 5„^m = (1 — f^)~ ^(u;,t;) into Eq. ( ^ ) and obtain the following 
equation for fuinction ^(c(;,f) 

(l-^;^)a^;^; -2{m + l)vg^ + {A™ - c^ t;^ - m(m + 1)} ^ = (4.6) 

Prom Eq. (^^) and Eq. ( [4.5D the set of the eigenfunctions for the generating function Z„ is equal to 

Z„ = (l-t;2)l 5„,i(i^,i;) et^i-^^) (4.7) 

The normalization of function Sn^i is given by the following equations |^] 

\\Sn,i{c,v)f = t dv \Sn,,{c,v)\'' = ^^llili) (4.8) 



J -I 2n + l 

The Green function of Eq. ( [4.4D has the form 
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(4.^ 



It should be stressed that the Y dependence of the Green function is determined by the condition 
En — Emin = —^n + ^0 = ^n ™ accordance with the Schroedinger equation (see Eq. ( 2.38| )). Eq = —loq 



is the lowest energy level which corresponds to the asymptotic state of Eq. ( |2.43| ). Therefore, Eq. (|4.9D 
explicitly demonstrate that the generating function Z tends to Z ^ Zasymp at high energies. 

This Green function leads to the following formula for the generating function Z{Y, u) 

Z{Y,u) = u + (4.10) 

u{l-u) ^^ ^ (l - e-^'r.Y\ 5„,i(if,l-2n) /^, ^^5^,1 (if, (- k ^T^ e" t (!"«) 



VI - (1 - 2ur ^^\ K J ||5„,i(i|,l-2^)P 



Fortunately, we have a rather detailed knowledge of the spheroidal functions in mathematics including 
the spectrum of A^ (see Refs. |2C, ^]). The series in Eq. ( 4.10| ) is well converged and, practically, we need 



to take into account only about 10 terms. 

Eq. ( 4.10| ) can be easily generalized for arbitrary initial condition: Z{Y = 0;u) = F{u). 



Indeed, the solution has the form 



Z{Y,u) = F{u) + — ^iL^Lx (4.11) 

^ ' ^ ' ^1 - (1 - 2uf ^ ' 



^^"E 



e-^ 



IY\ Sn,i{i^^,l-2u) f\ dCSn,i{q,C) [VT^L(^F{^)')e^ 



f (1-0 



"""' ^iV ^- J ||5„,i(if,l-2«)|P 

where L is the Sturm - Liouville operator defined in Eq. ( |2.30| ) . In coming papers we will discuss the master 
equation with different initial conditions that reflect the interaction with nuclei and other targets E3, H5[ 



5. Numerical estimates 

In our numerical calculations we use the package for mathematica written by P. Falloon (see Ref. 
which allows us to compute both the eigenvalues A^ and the spheroidal wave functions Sn,i with fixed c 
given by Eq. (|4.5|). The result of our calculations at fixed k = 16 which is related to the QCD coupling 
constant Ug = 0.25, are plotted in Fig. ^ and Fig. |8[ 

Fig. 1^ shows the generating function as function of u at different energies. One can see that at 3^ ^ 1 
the solution approaches the asymptotic solution of Eq. ( 2.43| ) ( the low dashed line in Fig. |^ in accordance 



with our expectations [42, 41, 13]. To understand the scale of the energy dependence we need to go back 
to the rapidity Y from y. As we have discussed 

Y = ^,/ ^, y = -^ — y (5.1) 

r(l ^ 2) ai.UJBFKL 
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Figure 7: The u dependence for the generating function Z{Y^u) (Fig. M-a) and the scattering amplitude N{Y,u) = 
1 — Z{Y,u) (Fig. u[h) at different values of rapidity. In this figure we plot Y = y (see Eq. (2.2) ). At y = oo the 
asymptotic solution of Eq. (2.43) is shown. k = 16. 



where ujbfkl is the intercept of the BFKL PomeronQ, namely, ojbfkl = 4 In 2 05. For as = 0.25 
which we use in this paper, Y = 23 y. Therefore, y = 0.3 gives Y — Yq = 6.9 (Iq is the initial rapidity 
from which we can apply our high energy approximation). 

Fig. ^ shows the main result of this paper: the energy (rapidity) dependence of the scattering ampli- 
tude. One can see that the exact solution approaches the asymptotic solution as has been foreseen in the 



approximate solutions to the problem (see, for example, Refs. [42, 41, |T^). 

This result eliminates the last uncertainty that the high energy limit of the scattering amplitude does 
not show the black disc behaviour but we have rather gray disc one as has been discussed before in Ref. 
[ |13| |. It turns out that at 3^ > 0.5 (Y > 12.5) the difference between the exact solution and the asymptotic 
one is very small. However, for y < 0.5 the difference between the semi-classical approach and the exact 
solution is rather large and all approximate methods should be taken with the justified suspicion in this 
kinematic region. 

In Appendix D we develop the numerical method specially for small values of Y and we compare it 
with the exact solution of Eq. (|4.10| ) . 



6. Comparison with approximate approaches 



6.1 Exact solution and the semi-classical approach 

As one can see from Fig. |^ the semi-classical approach preserves the singularity structure of the exact 
solution and correctly reproduces the high energy limit of the amplitude at any value of u but it cannot 
describe the approaching to this limit. The reason for such a behaviour is clear if we compare the spectrum 
A„ at low n of the semi-classical approach given by Eq. ( p. 22 ) with the values of A„ of the exact solution 
(see Table 1). In the semi-classical solution at small values of n AA„ = A^+i — A„ = k while in the exact 
spectrum AA„ < k. 
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Figure 8: The rapidity dependence of the scattering amplitude N{Y, u) at different values of u. The dashed curves 
are the semi-classical solution of Eq. ( |3.26 ), while the solid lines are the exact solution to the master equation. In 
this figure we plot Y = y (see Eq. (2.2) ) and k = 16. 



6.2 Exact solution and the mean field approximation (MFA) 



For the mean field approximation we need to neglect in Eq. {2A) the second term in the r.h.s. of the 
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equatterniThe solution is very simple (see Refs.|^, |lO|), namely 

-kY 



Z{Y-u) 



ue 



I + u {e-'^Y - 1) 
In Fig. y we plot the exact solution (solid lines) and the mean field approximation (see Eq. (pj])). 
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Figure 9: The rapidity dependence of the scattering amphtude N{Y,u) at different values of u. The dashed curves 



are the mean field approximation given by of Eq. (6.1), while the solid lines are the exact solution to the master 
equation. In this figure we plot Y = y (see Eq. (2.2) ) and k — 16. 



One can see that the mean field approximation cannot describe the scattering amplitude for values of 
u close to 1, but it gives a good approximation in the region of n < 0.6 in the wide energy range. However, 
even for large values of u the MFA fails to reproduce a correct approaching the scattering amplitude to 
the asymptotic limit. 

The physical reason for such a behaviour is clear. The typical value for the scattering amplitude 
at low energy 7 = 1 — n ~ 1/k < 1 for the dipole-dipole scattering. Therefore for the dipole-dipole 
scattering the exact solution leads to a different behaviour than the mean field approximation. However, if 
we have the interaction of dipole but with the nucleus target the amplitude at low energies is proportional 
to As 'y (dipole — dipole) = As / n = 1 — n ps 1 (see for example the paper of Kovchegov in Ref. ||T7|). 
The second case when the mean field approximation can work is the deep inelastic scattering. Indeed, we 
can model this process assuming that k = l/a|((5g). However, the scattering amplitude for low energies 
is proportional to a|.(l/i?^) where R is the target size. Therefore, for the running QCD coupling we can 
assume that 1 — u = 7 = a|(l/i?^) 
deep inelastic processes |]l|]. 



1 > 1/k. and the mean field approximation can describe the 
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6.3 Exact solution and the Mueller-Patel-Salam -lancu (MPSI) approximation. 



It was shown by Mueller and Patelp^] ( see also Refs. |37, |2^, |3^, ^) that in a restricted region of energies 
we can consider the Pomeron loops larger than y/2. Such loops we can sum and the answer for our model 
in ouir notation looks as follows (see Ref. p6|) 



N{Y, u) 



'l-u)e 



kY 



exp 



'1 — u) e 



kY 



r 0, 



(1 — u) e 



kY 



(6.2) 



Fig. ^ shows the comparison of the solution given by Eq. ( |6.2| ) (dotted lines) with the exact solution (solid 
lines). One can see that only in the very limited range of energies Y < 0.1 and at small values of u the 
MPSI approximation can describe the scattering amplitude. 
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Figure 10: The rapidity dependence of the scattering amplitude N(Y,u) at different values of u. The dashed 
curves are the MPSI approximation given by of Eq. (|6.2[), while the solid lines are the exact solution to the master 



equation. In this figure we plot Y = y (see Eq. (2.2)) and k = 16 
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6.4 Exact solution and the gray disc behaviour. 

One of the most interesting results of this paper is the asymptotic gray disc behaviour at uhra high energies: 
the solution does not go to 1 but to a constant which is smaller than 1 and which depends on the initial 
conditions or, in other words, on the value of the scattering amplitude at low energies. Such a solution 
does not contradict to the unitarity constraints and it is consistent with some of the previous approaches 



see Refs. [42, |l^). However, among the approximate approaches to the solution of the master equation 



situations is quite different: most of them: the mean field approximation, the MPSI solution as well as the 



solution of Ref. |15] lead to the black disc behaviour and only the semi-classical approach gives a constant 
limit at high energy that is less than 1. Therefore it is worthwhile to discuss the gray disc limit in more 
details. 

The scattering amplitude is equal to [^, 11] 



n\ 



NiYMYo)) = 1 - Z(y, 1-7(^0)) = E(-^)'^'('^o)[Y.Pn{Y) 

fc=l \n=l '^ ^' , 



7(Yo)<l 

n=l 

where ■j{Yq) is the scattering amplitude of one dipole with the target at low energy. If this amplitude is 



7(>o) 5] nP„(y) = (|n|)x7(yo) (6.3) 



very small Eq. ( |6.3| ) shows that the scattering amplitude for the dipole - target scattering at high energy 
is equal to the product of the low enegy amplitude by the average nuenergy((|n|)) of dipoles at low energy. 
Therefore, to obtain the amplitude that equals to 1 at high energy we need to assume that (|n|) becomes 
large (at least of the order of 1/7(1^0) )• Large number of low energy dipoles is natural to expect in the mean 
field approximation where the tree ('fan') diagrams contribute (see Fig. ||). On the other hand the mean 
field approximation violates the i-channel unitarity or, in other words, it gives an incorrect normalization 
of the parton wave function (see Ref. 1^31 )• Therefore, in the correct answer we rather expect that the 
average number of partions should be of the order of 1. Indeed, the enchanced diagrams of Fig. |l| are 
the most essential at high energies to protect the t - channel unitarity |^. In these diagrams we expect 
the average number of dipoles at low energy to be the same as in our initial condition, namely, (|n|) ~ 1. 
In this situation the behaviour of the amplitude at high energies depends on the value of the scattering 
amplitude at low energies as we see in our exact solution. 

This result we could also foresee if we considered the first corrections to the vertex of splitting of one 



Pomeron to two Pomerons (transition P — > 2P, see Fig. 11). First thing that we can see from this picture. 



which graphically shows the result of integration with respect to y2, that the first corrections lead to an 
increase of the value of the 2P -^ P vertices with rapidities yi ( the first diagram in Fig. 11) and y ( the 
second diagram in Fig. ^). Indeed, the simple diagram in Fig. |ll| has the following explicit expression 

A{Fig.^ = r2(1^2)r(2^1) / dyi T dy2 e^(^^2){y+?/i+?/-2/2) 

Jy Jo 

= r(i^2)r(2^i) r dyi Lni-^^){Y+y.) _ ^ni-.2) (Y+y,+y)\ ^54) 

J V 



25 




Y=0 
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Y=0 



Figure 11: The first correction to the P -^ 2P vertex due to the merging of two Pomerons into one Pomeron. Solid 
lines are Pomerons. Lines with arrows denote the interacting dipoles. 



Having this observation in mind we solve the mean field equations but assuming that the vertex for 
P ^ P transition (A) is not equal to the vertex of splitting of one Pomeron to two Pomerons r(l -^ 2). 
The equation has the form 



dZ{Y;u) 



-An + r(l 



2)u' 



dZ{Y;u) 



dY "" ' ^ ' ' du 

Eq. (|6.5D has an obvious solution, which for the scattering amplitude has the form 



N{Y,^{Yo)) 



7(lo) e^ ^""-^"^ 



Y-Yo > 1 



A 



1 + 7(yo)r(i^2) {e^o^-yo) - i) 



r(i ^ 2) 



(6.5) 



(6.6) 



One can see that if r(l ^^ 2) > A we expect a gray disc. It is easy to check that the MPSI approximation 
with different vertices for P ^ P and P -^ 2P transitions give the same answer as Eq. ( |6.6| ). Therefore, 
Eq. ( |6.4D shows that the gray disc behaviour is very likely for the exact solution. 

Formally speaking we expect to have two asymptotic solutions (Zoo) which stem from the master equa- 
tion with dZoo/dY=0, since the equation is the second order differential equation. The linear combination 
which satisfy the boundary conditions Z^o^u = 0) = and Zoo{u = !) = !, is the solution to the problem 
at high energy, which leads to a gray disc behaviour of the scattering amplitude. These two boundary 
conditions are very natural since they say: (i) for u = , that if amplitude is equal to 1 at y = none 
of interactions can change this fact due to the s-channel unitarity; and (ii) for u = 1, that starting with 
zero amplitude for dipole interaction at y = 0, we cannot obtain an amplitude different from zero in an 
interacting system. Therefore, neither the fact of two asymptotic solutions nor the boundary conditions 
depend on the particular form of interaction. For example, they do not depend on whether we took into 
account the possibility for two Pomerons to re-scatter in two Pomerons as it is done in this paper, or 
whether we neglect such interaction as was assumed in Ref. [^. In both cases we have the asymptotic 
solution, that satisfies the boundary conditionssatisfiesform of the gray disc. One of the main result of 
this paper is the proof that nothing unusual happens with the system in the process of evolution from 
y = to very large values of Y and the resulting solution approaches the asymptotic one. To understand 
this result it is instructive to share with our reader what kind of an alternative solution we could expect. 
Frankly speaking, we can expect in such kind of problem the asymptotic solution that has a singularity 
in the singular points of this equation (u = and u = 1). Even more, the solution of the mean field 
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approximation, for which we have the same boundary condition, has an asymptotic solution, which we can 
write in the form 

Z^{u) = e{l-u) (6.7) 

One can see that Z^[u) = for n < 1 and Zoq{u) = 1 for n = 1. However, this solution has an singular 
behaviour and Z'^ = (5(1 — u) . Nevertheless, since (1 — u) (5(1 — n) = this is a solution. The solution 
that we found shows that in general case there is no such singular solution. Of course, the last claim is 
based on our believe that the physical problem has one and the only one solution. One can check that 
(1 - u){d^Q{l - u)/(f'u) ^ zoo at n ^ 1. 

Our conclusion is very simple: the asymptotic behaviour of the exact solution is not only the gray disc 
but this kind of behaviour gives us a check how well an approximate method is able to discuss the high 
energy scattering. However, it should be very carefully stated that the real QCD case can be different; 



and we have arguments for both scenarios: black disc behaviour |48, 49 1 and the gray disc one [13|. The 
difference is in the extra integration over sizes of produced dipoles. The scattering amplitude in the Born 
approximation of perturbative QCD is very large for the dipoles of the large sizes. Therefore, the equation 
of Eq. ( |6.3D - type can give the black disc even if the average number of dipoles is about 1 since it can 
be one dipole but of very large size. However, the question, whether we can trust the perturbative QCD 
calculations for such dipoles , is the open question. 

7. Conclusions 

The core of this paper is the exact analytical solution to the BFKL Pomeron calculus in the zero transverse 
dimension given in section 4. The method, that we proposed, is very general and can be applied to a 
various kind of practical problems in future. 

The scattering amplitude at high energies, which stems from the exact solution, shows the gray disc be- 
haviour of the asymptotic solution (see Eq. ( |2.43| )). Therefore, we filled the hole in the proof demonstrating 
that the asymptotic behaviour of Eq. ( 2. 431) satisfies the initial conditions. 



We found thainitialapproximations to the master equation, developed in the past (see Refs. |^, ^, 



13, 14, [T3I), work quite well at large values of rapidity iY > 10) and can be used in either asymptotic 
estimates or for an illustration of the main properties of the high energy interaction. The advantage of all 
of them is their simplicity and transparent physics ideas that are behind these methods. In this respect 
the semi-classical approach is especially attractive since it gives the clear understanding of the spectrum 
of the problem and the behaviour of the eigenfunctions. However, the comparison with the exact solution 
shows that we cannot trust the approximate methods for low energies, namely, for Y < 10 where there 
exist the most of the experimental data including ones that we anticipate from the LHC. 

We firmly believe that the solution for the BFKL Pomeron calculus at zero transverse dimension is the 
necessary step in our understanding of the behaviour of the scattering amplitude at high energies in QCD. 
The real QCD problem is much more complicated but the fact, that we can compare the approximation 
methods with the exact solution for our toy model and found how they work, might open a way to solve 
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the QCD problem at high energies developing different approximations (see Refs.[^, 29, 3C, ^] where one 
can see the first attempts of such kind of approaches) . 
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A. Appendix A: General way for solving linear nonhomogeneous boundary problem. 



At this appendix we describe the general way for solving linear non-homogeneous boundary value problem 
using the appropriate Green's functions, that can be constructed from the eigenfunctions of the Sturm- 
Liouville problem. 

Here we are dealing only with the issues related to solutions of the master equation, which is equivalent 
to the BFKL Pomeron calculus in zero transverse dimension, while for more detail information, related to 
this subject, one can find in Ref. [ p2| . 

First of all, it is easily to check this using trivial transformations: 



s{v) 



f{v] 



exp 



q{v) 



ajv) 

fiv) 



dv 



fiv) 



exp 



p{v) = exp 
9{v^ 



ajv) 

fiv) 



dv 



fiv) 



dv 



the general linear non-homogeneous partial differential equation of the parabolic type: 



dy 



d'^Z dZ 



(A.l) 



(A.2) 



can be reduced to the equation which is suitable finding a solution to the Sturm-Liouville problem: 



s{v) 



dz_ 

dy 



d_ 
dv 



p{v) 



dz_ 

dv 



{v)Z + s{v)ip{y,v) 



(A.3) 



Denoting by ^{y,v) = s{v)(p{y,v) we have the equation: 



dZ_ _ d_ 
dy dv 



■ dZ' 



q{v)Z + <^iy,v) 



(A.4) 
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For arbitrary initial Z{0,v) = 6{v) and boundary Z{y,vi) = a, Z{y,V2) = P conditions the general 
solution is given by the following expression (see Refs. |2^, p2|): 

Z{y,v) = / ^T,0-G{y-T,v,OdTd^+ s{0 ■ 9(0 ■ G{y,v,0 d^ 

Jo J vi Jvi 

ry ry 

+ p{vi) a- Ai{y-T,v)dT + p{v2) P ■ A2{y - t,v) dr (A.5) 

Jo Jo 

Here the Green's function is given by: 

Giy,v,0 = f;^4|^7%^exp(-A„y); ||Z„||2 = H s{v) ■ Zl{v) dv ■ (A.6) 

„=1 ll^nll Jv^ 

where the A^ and Zn{v) are the eigenvalues and corresponding eigenfunctions of the following Sturm- 
Liouville problem for the second order linear differential equation: 

{p{v)Zy+[Xsiv)-q{v)]Z = (A.7) 

with the boundary conditions 

Z{vi) = a, Z{v2)=(3. (A.8) 

The functions Ai A2 are expressed via the Green's function: 

Ai{y,v) = lim —G{y,v,0; Miy,v) = lim —G{y,v,0- (A.9) 

(^vi at, ^^V2 ot 



The most complicated part in Eq. ( [A.5| ) is the two last terms that include functions Ai, A2 and require 
applying of some regularization technique in calculating them. Fortunately, we can find a simple transfor- 
mation after which these two terms will disappear. Namely, subtracting from Z(y, v) the initial conditions: 

Ziy,v) = Z{y,v)-Z{0,v) (A.IO) 

we will obviously have equation for function Z{y, v) with zero boundary conditions: 

Z{y,v^) = Z{y,v^) = (A.ll) 

and zero initial conditions: 

Z{0,v)=0. (A.12) 

The only price for this transformation is transformation of the free term in Eq. ( |A.4| ): 



d 

Hy,v) = '^{y,v) + — 



■ .dz{o,v) 

p{v) 



dv 



q{v)Z{{),v) (A.13) 
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Thus we have only the first term (term which has nonzero contribution) in Eq. ( [A.5| ) for the function 
Z{y,v): 



Z{y,v) 



y [■'"2 



<l>iT,0-Giy-T,v,^)drdC 



(A.14) 



'0 Jvi 

and final answer obviously has the following form: 

Z{y,v)=Z{0,v) + Z{y,v)=Z{0,v)+ / ^t,0 ■ G{y - T,v,OdT d^ 

Jo Jvi 



(A.15) 



Note, that using Eq. ( A.15| ) we reduce Eq. (A. 3) to the equation for Z{y,v) which has the same form as 
Eq. (A. 4) but with the inhomogeneous term $(y,t>) given by Eq. ( |A.13| ). The advantage of this transform 
is clear if we compare the rather compact answer of Eq. ( A.15| ) with the solution in the form of Eq. ( |A.5| ). 



B. Appendix B: The exact semi-classical solution. 

At this appendix we calculate exactly functions (p {u!,u) for the semi-classical solution (see section 3.1) 
and compare them with the approximations at small and large values of w, which were used for calculation 
in this section. 



Recall that 



2 I ^ ^ V^ ^ k2 n (1 - u) 



du 



Therefore, functions (j)^ (see Eq. (W^ and Eq. (p^)) are defined by the following integrals: 



(B.l) 



<^±(n) 



where the first integral is trivial: 



K 

- < 1 ± 

2 I 



'1 + 



4 u 



K? u' (1 — U') 



du = Ii + I2 



(B.2) 



Ii = - I ^du' = f(n-l) (B.3) 

The second term ^ is a bit more complicated and can be expressed using well known elliptic integrals: 
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K 

= ^2 


1: /' ^ 


4 u) 

k2 u' (1 - 


-— du' 
-u') 




/I 4u; 
V4+.2 


1 




1 + IQlo/k? 



EE ArcSin[l - 2u] 



1 + IQuj/k? 



**We very grateful to our referee, who noticed that I2 can be reduced to this relatively simple form. 



(B.4) 
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Figure 12: This plot shows the exact function 0+(a;,M) (thick dotted hnc) and the approximate function (^(w) of 

— 16 for the plot. 



Eq. (|3.6[) (thin solid line) for small values of uj (i.e. a; < k 



Where EE{(p,m) is the elliptic integral of the second kind |25] 



EE{ip,m) = / (1 
Jo 



m ■ sm 



1/2 



'i6)y''Ue 



(B.5) 



We can see that, in spite the fact that (p can be calculated explicitly (in radicals), the received 
expression is to complicated in order to use it for the calculations of the functions $ (w) in Eq. 



Nevertheless, this exact calculation very useful by two reasons: (i) using exact form of cj) we can estimate 
the accuracy of our approximate expression for functions </)^(w,n), and (ii) we can find the separation 
parameter between small and large co that we used in our semi-classical solution of section 3.2. 



The plots in Fig. 12 and in Fig. 13 illustrate how good our approach for small and large u>, which we 



used for calculations of - 



C. Appendix C: Hermitian Hamultotian description of rapidity evolution. 



At this appendix we give the answer to the important question: whether rapidity evolution in zero transverse 
dimensions can be described by hermitian Hamiltonian ? Or, in other words: does we have probability 
conservation in these processes ? 

First of all, let us rewrite our master equation Eq. (p.l[) : 






du dv? 



(C.l) 
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Figure 13: This plot shows the exact function 0+(a;,M) (thick dotted Unc) and the approximate function (^(w) of 
Eq. (|3.6|) (thin solid line) for large values of a; (i.e. w > k^). k — Vo for the plot. 



in the form of the Fokker-Plank equation, using the simple transformation: Z{y,u) = n(l — u) ■ Z{y,u) 



dZ d f ~ \ d^ 



It is natural to define the Fokker-Plank operator as 

Up ^ - ^a{u) + ^b{u) 
where in our case a{u) = ku{1 — u) and b{u) = n(l — u). Introducing potential <I>(u) such that: 



(C.2) 



(C.3) 



$(n) = ln[6(n)] 



we can rewrite Lpp in the following form: 



"a(w' 



ui b{u' 



— du' = ln[u(l — u)] — K 



d 



Lpp = -b{u)e 



du 



(C.4) 



(C.5) 



Now it is obviously that Lpp is non-hermitian operator. Nevertheless, operator Lh = e^^'^>''^Lppe~^^^'i'^ 
is herniitian (in some sense we restored the hermitian property of our system). Indeed we can check this 
by the direct calculation: 



(t)i{u)e^^''^LFp(t)2{u)du 



ui 



0,(n)e*(")^6(n)e-*(")^e*(")(A2(u)du 

Ml 



du 



du 
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A*,(«)e»'"' 



h{u)e 



-$(m) 



lHu)e-^" 



Ul 

ou ou 



du 



Ml 
W2 



(/)2(n)e*(")Li7p<^i(u)dn (C.6) 



Ml 



where <f)i{u) and 02 (^) are two functions, which satisfy boundary conditions. Therefore, we have demon- 
strated that fe*(")Li.p') = L],pe*(") = e*(")LFP 



It is important to note that transformation Lh = e '"" ^Lppe"^"" ^ can be ahernatively considered as 
a transformation of the appropriate eigenfunctions. If (/>„(u) are the eigenfunctions of the Fokker-Planck 
operator Lpp with the eigenvalues An, the V'n(^) = e*^""^ (/>„(«) wih be the eigenfunctions of operator Lh 
with the same eigenvalues. Actually, appropriate transformation of eigenfunctions was used in our paper 
in order to restore Hermitian property of our system. 

It is obviously that both of system 0„(u) and ipn{u) predefine the orthogonal and complete systems: 

• Orthogonality: 

• Completeness: 

e*^"^ En <t>n{u) ■ Mn') = En M^) ' ^n(n') = 5{u - u') 

D. Appendix D: The exact solution as a power series in rapidity. 



Our master equation Eq. (2.1) can be solved exactly for the arbitrary initial conditions: Z{Y = 0; u) = f{u) 
in terms of the power series of rapidity. This method has a number of exclusive advantages: (I) it applicable 
for all reasonable initial conditions, which allow us to apply such solution for the different scattering 
processes at high energies; (II) with this method we can easily solve all modifications of our master equation, 
which incorporates next order corrections with the only one restriction that the variable coefficients in the 
diffusion equation can be represented as power series in u); (III) this method extremely good for small 
values of rapidity y and in vicinity of u = and u = 1 (we will demonstrate this below with numerical 
calculations). The last argument especially important, since all solutions, which are based on approximate 
methods, lead to large errors just for small values of rapidity y and/or in vicinity of u = and n = 1. 

Unfortunately, solution in terms of power series of rapidity is converged relatively slowly, thus this 
method is not applicable for calculations of the asymptotic behaviors of the scattering amplitude. Never- 
theless, we believe that proposed method will be very useful for practical calculations, because most of the 
experimental data, available now, are concentrated at small or intermediate values of rapidity. 

It is very instructive to start with the direct decompositions of function Z[y,u): 
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z{y,u) = Y,Cn{u)-y 



fD.ll 



n.=0 



It is obviously the first term in such expansion corresponds to appropriate initial conditions, indeed at 
y = Q we have: 



z(3^ = o,u) = Y,Cn[u)-{y = dT = c, 



oW 



f{u) 



(D.2) 



n=0 



Now by substitution of Eq. (|D.1| ) into our master equation Eq. ( |2.lD and after arranging the appropriate 
powers of y we will get the following recurrence equations for the functions Cn{u): 



C„+i(u) 



n(l 



u) 



n + l 



d^Cn{u) dCn{u) 



du^ 



du 



(D.3) 



Therefore, once we know initial conditions of our problem Eq. ( D.2 ), all series expansions Eq. ( D.l| ) are 
well defined and can be calculated up to arbitrary order. Finally, introducing the differential operator L" 
as follows 



V[f{u)] = L V-'[f{u)] 



the answer can be written in the very elegant form: 



L 



u{\ — u) 
n + l 



dv? 



d_ 

du 



(D.4) 



z{y,u) = f{u) + Y.L''[f{y)]-y'' 



(D.5) 



n=l 



This result is nice from the analytical point of view: (1) it haven't singularities, (2) the expansion is well 
defined for all orders. Nevertheless, for numerical estimations probably not so usefully. The reason for this 
is clear: series expansion Eq. ( p.5| ) is converged relatively slowly and even more, the radius of convergence 
is the function of u. For example, taking first 20 terms of our expansion we can ensure solution Z(Y,u) in 
rapidity rang Y G [0,5]. So from numerical point of view this solution can be considered as solution for 
small and intermediate values of rapidity. 

Fortunately, we can improve numerical property of proposed solution by simple modification, namely 
we are going expand not the function Z{y, u), but function <p{y, u) = ln[Z{y, u)]. It is obviously that this 
modification was inspired by our experience in semi-classical solution, but now we are going to keep all 
terms and propose exact solution. 

The first step, which we do is rewriting our master equation Eq. ( |2.l[) in terms of function ip{y, u) = 
ln[Z{y,u)], by simple substitution we get the following equation for Lp{y,u): 
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Now one again, we search a solution for this equation in terms of power series of rapidity: 

oo 

if{y,u) = ^Bn{u)-y^ (D.7) 

n=0 



Substituting Eq. ( p.TD into Eq. ( D.6| ) and arranging the appropriate powers of 3^ , we obtain the following 



recurrent relations for functions Bn{u): 



R t \ u{l-u) 



dBk{u) dBn-k{u) , d^Bn{u) dB^iu 



r)ii r)ii 



K- 



du du du^ du 

k=0 



n + 1 
The appropriate initial condition for function (p{y, u) is also simple In transformation: 



(D. 



ip{y = 0,u) = J2Bn{u)-{y = 0r = Bo{u) = ln[/(n)] (D.9) 

n=0 

The last solution numerically is much more stable and, taking first 20 terms of our expansion, we 
can ensure solution Z(Y,u) already in rapidity rang Y G [0,8]. It is important to note that proposed 
logarithmic transformation works well in all cases, except the specific initial conditions: Z{y = 0, u) = 
const (fortunately, for our physical environment such initial conditions does not look reasonable). 

Finally, we demonstrate the quality of this method comparing it with the exact solution of Eq. ( [4.10 ) . 
It is obviously that for small values of rapidity this method work quite well. 

E. Appendix E: The interaction procedure for the spectrum of the problem. 

In this appendix we show that the iteration procedure that has been described in 2.3, works quite well 
at least for calculations of the eiegenvalues for the master equation (see Eq. ( |2.1| )). Starting with the 
exact spectrum of the equation at large values of n (see Eq. ( p.l7| )) and using the set of eigenfunctions we 
can develop the regular perturbative approach based on Eq. ( p.44| ) with Uq (G) being the rectangular-well 
potential and choosing V (0) = ^ sin^ 6. The corrections to the eigenvalues (AA„) are determined by the 
simple formula 



AAn = / dG V{Q) \^{Q)Y (E.l) 

JO 

and the resulting 

A„ = n2 + AA„ (E.2) 



35 




0.2 0.4 0.6 0.8 
U 



Figure 14: In this plot we compare the exact solution of Eq. ( 4.1C ) - a solid line and the solution, given by the power 
series in rapidity, taking into account first 20 terms (see Eq. ( D.5 )) - a dashied line. This comparison is performed 
for two different values of rapidity. As we can see from the plots these solutions are close. In principle, the exact 
solution is a bit steeper then the power series solution, which obviously reflects the fact that we took into account a 
limited number of terms in this expansion. 



Eigenvalues 


Ai 


A2 


A3 


A4 


A5 


Ae 


A7 


As 


A9 


Aio 


Spheroidal (exact) 


13.0 


14.13 


21.40 


28.63 


38.45 


50.32 


64.24 


80.19 


98.15 


118.13 


Asymptotic (A^ = n ) 


1 


4 


9 


16 


25 


36 


49 


64 


81 


100 


First iteration 


14.25 


14.75 


21.25 


29.75 


40.25 


52.75 


67.25 


83.75 


102.25 


122.75 



Table 1: The comparison of the first ten eigenvalues of the master equation with the asymptotic estimates and with 
the first iteraction described in section 2.3. 



One can see from the Table I that this first iteration leads to the eigenvalues which are very close to the 
exact eigenvalues of the master equation. 

Prom this table we see that the first iteration reproduces correctly the values of the eigenvalues at 
small n. Therefore, we can conclude that the iteration procedure, described in section 2.3, gives the way 
to approach the exact solution with a good accuracy. 
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